Draft version September 13,2012 

Preprint typeset using style emulatcapj v. 12/16/1 1 


THE RISE AND FALL OF STAR FORMATION HISTORIES OF BLUE GALAXIES AT REDSHIFTS 0.2 < z < 1.4 

Camilla Pacific!, 1,2 Susan A. Kassin 34 Benjamin Weiner, 5 Stephane Charlot, 1 2 Jonathan p. Gardner, 3 4 5 et al. 

Draft version September 13, 2012 

ABSTRACT 

Popular cosmological scenarios predict that galaxies form hierarchically from the merger of many progen- 
itors, each with their own unique star formation history (SFH). We use the approach recently developed by 
Pacific! et al. (2012) to constrain the SFHs of 4517 blue (presumably star-forming) galaxies with spectro- 
scopic redshifts in the range 0.2 < z < 1.4 from the All- Wavelength Extended Groth Strip International 
Survey (AEGIS). This consists in the Bayesian analysis of the observed galaxy spectral energy distributions 
with a comprehensive library of synthetic spectra assembled using state-of-the-art models of star formation 
and chemical enrichment histories, stellar population synthesis, nebular emission and attenuation by dust. We 
constrain the SFH of each galaxy in our sample by comparing the observed fluxes in the B, R, I and K 3 bands 
and rest-frame optical emission-line luminosities with those of one million model spectral energy distributions. 

We explore the dependence of the resulting SFH on galaxy stellar mass and redshift. We find that the average 
SFHs of high-mass galaxies rise and fall in a roughly symmetric bell-shaped manner, while those of low-mass 
galaxies rise progressively in time, consistent with the typically stronger activity of star formation in low-mass 
compared to high-mass galaxies. For galaxies of all masses, the star formation activity rises more rapidly at 
high than at low redshift. These findings imply that the standard approximation of exponentially declining 
SFHs widely used to interpret observed galaxy spectral energy distributions is not appropriate to constrain the 
physical parameters of star-forming galaxies at intermediate redshifts. 

Subject headings: galaxies: evolution - galaxies: star formation - galaxies: stellar content 


1. INTRODUCTION 

Constraints on the stellar content of galaxies are often de- 
rived from multiwavelength observations by assuming that the 
star formation history (SFH) of an individual galaxy can be 
approximated by a simple declining exponential function of 
the form er t,T , where t is the galaxy age and r the star for- 
mation timescale. Such “r-models” have been used success- 
fully to estimate, for example, the stellar masses of nearby spi- 
ral galaxies from fits of rest-frame optical/near-infrared colors 
(e.g.. Bell & de Jong 2001). Despite this success, an increas- 
ing number of analyses have pointed out the limitations of 
this approach, particularly for applications to studies of high- 
redshift galaxies. This is the case of Lyman-break galaxies at 
redshifts z > 2 (Papovich et al. 2001) and blue star-forming 
galaxies at z ~ 2 (Shapley etal. 2005), for which the addi- 
tion of another stellar component on top of a r-model can 
change the derived stellar masses by a factor of several. In 
fact, the major episodes of star formation in ultraviolet (UV)- 
luminous Lyman-break galaxies could last only a few hun- 
dred million years (Stark et al. 2009). Reddy et al. 2012 also 
find that, for galaxies at 1.5 < z < 4, fits of broadband 
spectral energy distributions (SEDs) using exponentially de- 
clining SFHs produce star formation rates (SFRs) 5 to 10 
times lower than inferred from more direct estimates based on 
the combined analysis of the UV and mid-infrared emission. 
They conclude that SFHs in which the SFR rises in time are 
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more appropriate than declining r-models for these galaxies. 
This is reinforced by the conclusion by Maraston et al. 2010 
and Pforret al. 2012 that SED fits of galaxies at z ~ 2 us- 
ing declining r-models produce unrealistically low ages, and 
that exponentially rising r-models should be preferred for 
these galaxies. In fact, simulations of galaxy formation in 
a hierarchical universe predict complex SFHs, which are not 
well approximated by simple declining exponential laws (e.g., 
Lee et al. 2009; Wuyts et al. 2009). 

To better characterize the SFHs of galaxies at moderate 
and high redshift, we require more sophisticated, physically 
motivated spectral analysis tools. In this Letter, we achieve 
this goal by appealing to the approach recently proposed 
by Pacifici et al. 2012 to constrain the SFHs of 4517 galax- 
ies with secure spectroscopic redshifts in the range 0.2 < 
z < 1.4. Specifically, we fit the observed broadband SEDs 
and emission-line strengths of these galaxies using a com- 
prehensive library of synthetic spectra assembled by combing 
state-of-the-art models of stellar population synthesis, nebu- 
lar emission and attenuation by dust with star formation and 
chemical enrichment histories derived from the semi-analytic 
post-treatment of a large-scale cosmological simulation. This 
approach allows us to characterize the SFHs of galaxies in dif- 
ferent mass and redshift ranges in a statistically reliable man- 
ner. 

The paper is organized as follows. In §2, we present ob- 
servations of the B, R, I and K s fluxes and rest-frame op- 
tical emission-line luminosities of 4517 blue galaxies from 
the All-Wavelength Extended Groth Strip International Sur- 
vey (AEGIS, Davis etal. 2007). In §3, we describe the li- 
brary of synthetic SEDs built to interpret these observations. 
We report our results in §4 and present our conclusions in §5. 
Throughout this Letter, we adopt a standard ACDM cosmol- 
ogy with Gm = 0.3, Qa = 0.7, h = 0.7, and a Chabrier 
initial mass function (Chabrier 2003). 
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2. DATA 

We require our observational sample to cover a large range 
in galaxy stellar mass and SFR. In addition, to derive good 
constraints on galaxy SFHs and stellar masses, we need 
solid spectroscopic redshifts and accurate photometry span- 
ning the 4000A break to the rest-frame /-band. Further- 
more, emission-line luminosities are essential to accurately 
measure the current SFR (e.g., Brinchmann et al. 2004). The 
AEGIS Survey contains galaxies observed in a large wave- 
length range from X-ray to radio. From this catalog, we 
extract 6246 galaxies which have photometry at B, R, I 
(Coil etal. 2004) and K a (Bundy etal. 2006). The resulting 
sample is magnitude-limited at K a = 22.5. From this sample, 
we select 4517 potentially star-forming galaxies (U - B < 
1.0) that do not show contamination by an Active Galactic 
Nucleus (AGN). In particular, we reject the sources that are 
detected by Chandra and that are identified as AGNs fol- 
lowing the criteria by Kauffmann etal. 2003 (rest-frame op- 
tical emission-line ratios), Yan etal. 2011 (rest-frame opti- 
cal emission-line ratio versus color), and Donley et al. 2007 
(slope of the Spitzer IRAC SED). For all 4517 galaxies, 
Keck/DEIMOS spectra are available from the DEEP2 Red- 
shift Survey (Newman etal. 2012). They cover the wave- 
length range 6500-9100A at a resolution of 1 .4A full width 
at half-maximum (FWHM). The accuracy of the redshifts 
is 30 km s -1 . Such resolution allows us also to extract 
reliable emission-line fluxes of [O ll]AA3726.0, 3728.8, H/3, 
[6 IIIJA5007, and Ha when they fall into the observed wave- 
length range (Table 1). This set of photometric and spectro- 
scopic data covers rest-frame wavelength ranges from A ~ 
3300A-1.8/xm at z = 0.2 to A ~ 1650-9150A at z ~ 1.4. 

3. MODELING APPROACH 

We derive precise constraints on the SFHs and stellar 
masses of galaxies using a sophisticated tool based on phys- 
ically motivated models of SFH, stellar population synthesis, 
nebular emission and attenuation by dust. In this Section, we 
describe our library of models and how we fit the data. 

3.1. Library of model spectral energy distributions 

We build a comprehensive library of model star forma- 
tion and chemical enrichment histories by performing a 
semi-analytic post-treatment of the Millennium cosmolog- 
ical simulation (Springel etal. 2005) using the models of 
De Lucia & Blaizot 2007. Following Pacifici et al. 2012, we 
use a sample of 100,000 galaxy SFHs from z — 127 to the 
present time. We then enlarge the model library by draw- 
ing 10 different realizations of the redshift of observation 
(0.1 < z < 1.5) and of the evolutionary stage for each 
galaxy (see Sections 2.1 and 3.1.2 in Pacifici etal. 2012) to 
account for the fact that galaxies can be observed at differ- 
ent stages of evolution. For each realization, we also resam- 
ple the current (i.e„ averaged over a period of 10 Myr be- 
fore a galaxy is observed) specific SFR and gas-phase oxygen 
abundance, in the ranges —2 < log(ips/ Gyr* 1 ) < 1 and 
7 < 12 + log (O/H) < 9.4, respectively, to avoid biases in 
parameter retrieval. 

From this library of one million model SFHs, we generate 
SEDs by computing the emission from stars and gas. De- 
tails of how this is done are in Pacifici et al. 2012. The mod- 
els rely on the latest version of the Bruzual & Chariot 2003 
stellar population synthesis models and compute the nebu- 
lar emission using the photoionization code CLOUDY as in 


Chariot & Longhetti 2001 . We include the attenuation by dust 
with a two-component model as in Chariot & Fall 2000. For 
the dust model, we randomly choose the total optical depth of 
the dust (tv) between 0 and 3. Also for the dust model, we 
vary both the slope of the attenuation curve in the interstellar 
medium and the fraction of optical depth contributed by the 
dust in the interstellar medium in the ranges 0.4- 1.1 and 0.1- 
0.7, respectively. Finally, to simulate the absorption by the 
intergalactic medium as a function of redshift, we adopt the 
prescription of Madau 1995. 

This procedure allows us to go beyond the typical approach 
of using simple idealized functions at fixed metallicity to de- 
scribe the SFHs of galaxies. However, one caveat is that our 
results may still depend on the star-formation prescriptions 
adopted in the semi-analytic model. 

3.2. Fitting procedure 

We use the final library of 1 million model galaxies de- 
scribed above to measure physical parameters for the 4517 
galaxies with the multi-wavelength data presented in §2. We 
use a Bayesian approach, detailed in Pacifici etal. 2012. In 
brief, for each galaxy in our observational sample, we com- 
pute the likelihood that each model in the library describes 
the observations. Specifically, we fit to photometric fluxes in 
the observed B, R, I and K a bands, in addition to emission- 
line luminosities of [Oil], H/3, [O III] and/or Ha. The ex- 
act lines fit depend on the rest-frame wavelength coverage of 
the spectra, as given in Table 1 . Since galaxies in our ob- 
servational sample have accurate spectroscopic redshifts, we 
compute likelihoods only for the model galaxies built within 
0.05 of the spectroscopic redshift. We are thus fitting each 
observed galaxy with ~ 71,500 models. For each galaxy, 
we use the computed likelihoods to build probability density 
functions (likelihood distributions) for the stellar mass and the 
SFR normalized to the absolute observed A",, -band luminosity. 
The typical uncertainties (defined as half of the 16th-84th per- 
centile range of the cumulative probability density function) 
on the stellar masses and the SFRs are ~ 0.1 and 0.3 dex, re- 
spectively. For each galaxy in the sample, we also derive the 
best-fit SFH. 

4. STAR FORMATION HISTORIES OF BLUE GALAXIES 

Blue galaxies form a distinct sequence in the 
parameter space of stellar mass and SFR (see 
e.g., Brinchmann et al. 2004, Noeske et al. 2007, 
Daddi et al. 2007, Whitaker et al. 2012), the SFR increasing 
systematically with mass and redshift. This is referred to 
as the “star-formation main sequence” (Noeske etal. 2007). 
In this Section, we examine the average SFHs of low- and 
high-mass galaxies as a function of redshift. 

In Figure 1, we examine the best-fit SFH models for two 
example galaxies at similar redshifts, and which have simi- 
lar current SFRs (ip) but different stellar masses (M*). The 
first galaxy has a low stellar mass of M* = 4.3 x 10 9 M© 
and ip — 4.8M©yr -1 . The second galaxy is more mas- 
sive with M* = 2.9 x 10 lo M© and ip — S-SA^yr' 1 . It 
is interesting to note that the best-fit SFHs evolve inside the 
main sequence for their entire lifetimes. In the inset panels in 
Figure 1, the best-fit SFHs of these two galaxies are shown. 
These two SFHs show clearly the characteristic fluctuations 
over time that derive from the merging histories of the dark 
matter haloes and from the semi-analytic prescriptions used 
to model the baryonic component (gas infall, cooling, star for- 
mation, feedback). Neither of them would be well described 
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TABLE 1 

Median stellar masses and SFRs of AEGIS galaxies 


Redshift Observed emission lines log(M*/Af©) log[^>/(M©yr -1 )] number of galaxies 




Iow-mass“ 


low-mass' 1 

high-mass b 

low-mass 0 

high-mass 6 

0.20 < 2 < 0.45 

H/3,10 rri],Ha 

937 

10.48 

-0.185 

0375 

167 

16 

0.45 < z < 0.70 

H£,[Om] 

9.40 

10.56 

-0.065 

0.765 

272 

68 

0.70 < 2 < 0.85 

[Oll],H/3,[0 III] 

9.41 

10.54 

0.115 

0.855 

261 

124 

0.85 < 2 < 1.10 

[OII],H0 

9.43 

1052 

0.625 

1.155 

191 

105 

1.10 <2 < 1.40 

[Oil] 

9.50 

1054 

0.965 

1335 

96 

121 


“ Galaxies between 1 .6 and 4.0 X 10 a Af*/M© (blue dots in Figure 2) 

6 Galaxies between 2 -5 and 6.3 X 10 10 A/./M© (green dots in Figure 2) 



8 9 10 11 


log(H,/Uj 

Fig. 1 .— Measured SFR versus stellar mass for galaxies in the redshift 
range 0.85 < z < 1.1 are shown as small grey points. The blue dash-dotted 
line shows the best-fit model SFH of an example low-mass galaxy, while 
the green solid line shows the best-fit model SRI of an example high-mass 
galaxy. Inserted panels show the SFR versus lookback time for these two 
example galaxies. Both best-fit models have significant variations with time 
and neither would be well described by a r-model. 

by a r-model. 

In Figure 2, we assess the shape of the average SFH of 
galaxies as a function of stellar mass and redshift. In the col- 
umn on the left, we plot the star-formation main sequence 
for five redshift bins spanning 0.2- 1.4. We explore two stel- 
lar mass bins which are sampled at all redshifts: a low- 
mass bin (1. 6-4.0 xlO 9 M*/M 0 ) and a high-mass bin (2.5- 
6.3 xlO lc M*/M 0 ). Median stellar masses and SFRs in these 
two bins are given in Table 1 . In the middle and right-hand 
columns, we plot the average best-fit SFHs of the galaxies in 
the low- and the high-mass bins, respectively, and the root 
mean squares about these averages. To compute the average 
SFH, we place the best-fit model SFH of each galaxy at the 
redshift of observation. We then co-add the individual SFHs 
as a function of cosmic time and divide by the number of 
galaxies that form stars at that time. The small deviations 
in the SFHs of individual galaxies (as shown by the example 
galaxies in Figure 1) are generally washed out by taking an 
average of many SFHs. However, they are still visible in the 
top right panel (high stellar mass, low redshift) since there are 


only 16 galaxies in this bin. The number of galaxies in each 
individual plot is listed in Table 1 . 

From Figure 2 is clear that the average SFH of low-mass 
galaxies is characterized by a rising function at all redshifts 
with a progressively steeper slope with increasing redshift. 
This is likely because high-redshift galaxies have had less 
time to form stars than galaxies at lower redshift. We find that 
high-mass galaxies on average form more quickly than low- 
mass galaxies. They are characterized by bell-shaped. SFHs 
in which the SFR reaches a maximum and then declines. This 
is the case at all redshifts examined. From Figure 2 it is also 
clear that the lookback time at which galaxies have formed 
half of their total stellar mass shifts towards progressively 
higher lookback times with both increasing redshift of obser- 
vation and increasing stellar mass. For high-mass galaxies it 
coincides with the peak of their SFR. 

We make sure to include in our library a large variety of 
SFH shapes. These include shapes which are similar to de- 
clining, rising, bell-shaped and bursty. One way we have 
tested the robustness of our fits consists in examining the sec- 
ond and the third best-fit models for each galaxy. Doing so, 
the results shown in here do not change significantly. 

5 . DISCUSSION AND CONCLUSIONS 

Popular cosmological scenarios predict that galaxies form 
hierarchically from the merger of many progenitors, each with 
their own unique SFH. With the maturity of semi-analytic 
models of galaxy evolution, we now have the ability to model 
the hierarchical build-up of galaxies. Therefore, we no longer 
need to rely on simple analytic descriptions of galaxy SFHs, 
and can adopt a more sophisticated approach. 

In this work, we measure the SFHs of 4517 blue galax- 
ies in the AEGIS Survey using a sophisticated library of 
model galaxy SEDs. We build this library by performing 
a semi-analytic post-treatment of the Millennium Simulation 
(Springel et al. 2005) and including the latest stellar popula- 
tion models along with prescriptions for nebular emission and 
attenuation by dust. For each galaxy, we retrieve the best-fit 
model SFH and the probability density functions of its stellar 
mass and SFR. We reach the following conclusions: 

• The average SFH of high-mass galaxies rises and falls 
in a bell-shaped manner. This is inconsistent with sim- 
ple exponentially-declining or rising r-models. 

• Low- mass galaxies have rising and extended SFHs , also 
inconsistent with r-models. This implies that low-mass 
galaxies are currently more actively star-forming than 
high-mass galaxies. This is consistent with the trend in 
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Fig. 2.— Left-hand column: Star-formation main sequence (i.e., SFR versus stellar mass) for the observed galaxies divided into five redshift bins spanning 
02-1.4. Lew-mass galaxies (M* = 1. 6-4.0 x 10 9 Mq) are highlighted in blue, while high-mass galaxies (M* = 2.5-6.3 X 10 10 Mq) are highlighted in green. 
Middle column: Average best-fit SFH of the low-mass galaxies (black solid line) and root mean square about this average (blue shadow). Right-hand column: 
Average best-fit SFH of the high-mass galaxies (black solid line) and root mean square about this average (green shadow). The average SFHs of high-mass 
galaxies rise and fall in a bell-shaped manner. Low-mass galaxies have rising and extended SFHs, and are therefore currently more actively star-forming than 
high-mass galaxies. Both low- and high-mass galaxy average SFHs are inconsistent with simple idealized r-models. 


the star-formation main sequence (Noeske et al. 2007) 
and with the latest simulations by Behroozi et al. 2012. 

• The slopes of the average SFHs of both low- and high- 
mass galaxies are steeper at high redshift. This is likely 
because at high redshift, galaxies have less time to form 
stars compared with low redshift. 

• On average, the time at which half the stellar mass 
in a galaxy has formed (vertical black dashed lines in 
Figure 2) is a function of both redshift of observation 
and stellar mass. It shifts towards progressively higher 
lookback times with increasing redshift of observation 
and increasing stellar mass. This is the case for both 
low- and high-mass galaxies. 

• Both low- and high-mass galaxies evolve inside the 
star-formation main sequence for their entire lifetimes. 


Our conclusions rely on physically motivated models of 
SFHs that can be used to interpret any type of observations. 
The work presented in this letter will be extended by using ab- 
sorption lines in the spectra and observer-frame infrared pho- 
tometry to improve the constraints on the SFHs (Pacifici et al. 
in preparation). 
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